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Domain wall fermions and the strange quark mass * 

Matthew Wingate"" 

'^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 

The strange quark mass has been computed using a lattice action which possesses continuum-like chiral sym- 
metry to good precision, namely the domain wall fermion action. This talk surveys this action and the recent 
calculation of rria by the RIKEN/BNL/CU collaboration. This result is put into context by briefly summarizing 
other recent lattice studies. 

1. INTRODUCTION 

The strange quark mass is a fundamental 
parameter of the Standard Model, one whose 
value is poorly known. The two methods which 
have been used extensively to determine 
are QCD sum rules and lattice QCD simula- 
tions. Before one can compare results from 
these different methods fairly, a basic under- 
standing of each method's systematic uncertain- 
ties and approximations is necessary. This talk 
is focused on recent lattice calculations of 
emphasizing such uncertainties. The calcula- 
tion of TTis using domain wall fermions by the 
RIKEN/BNL/Columbia (RBC) lattice collabora- 
tion [0 is taken as a case study. The details of 
this work were presented at Lattice 99 and 
there was another description of the relevant lat- 
tice techniques here by Becirevic 1^ . 

2. DOMAIN WALL FERMIONS 

Symmetries are of fundamental importance in 
physics; in the case of QCD chiral symmetry 
is responsible for much of light hadron dynam- 
ics. However, chirally symmetric discretizations 
of the Dirac operator suffer from a "doubling" 
of the particle spectrum. For example replac- 
ing with a nearest -neighbor difference oper- 
ator leads to the free quark propagator S{p) — 
[z7^ sin(p^a) -I- m]^^ which, in the massless limit, 
has poles at all corners of the Brillouin zone 
G [0,7r]. Wilson added to the lattice action 
a second derivative-like operator which breaks 



Figure 1 . Schematic drawing of the light fermion 
wavefunctions in the extra dimension. For ap- 
propriate choice of the domain wall height M, 
a right-handed (RH) mode appears on the left 
boundary and a left-handed (LH) mode on the 
right boundary, TV^ sites away. The mixing be- 
tween the two modes within the extra dimension 
decays exponentially as TV^ increases. 
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chiral symmetry and gives the doubler states a 
mass inversely proportional to the lattice spacing. 
Another approach, by Kogut and Susskind (KS), 
is to stagger the spin degrees of freedom and in- 
terpret the doublers as different quark flavors. Al- 
though this method maintains a L/(l) remnant of 
the continuum chiral symmetry, flavor symmetry 
is badly broken. Both actions recover the sym- 
metries of continuum quarks only in the a 
limit. 

An exciting and fruitful suggestion for decou- 
pling the chiral and continuum limits started 
with Kaplan H] who pointed out to the lattice 
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community that a chiral 2A;-dimensional mode 
is bound to a mass defect (or domain wall) 
in a 2fc + 1-dimensional space. Specifically, in 
Shamir's boundary fermion variant j^,^, there is 
a massless right-handed (RH) mode stuck to the 
Ad boundary of a semi-infinite bd space, given an 
appropriate choice for the Sd mass M . When the 
fifth dimension is made finite, a left-handed (LH) 
mode appears on the other boundary (see Fig. |^). 
The wavefunctions of the light modes decay ex- 
ponentially within the fifth dimension, so there is 
a residual mixing between the RH and LH states 
which decreases as one increases Ng, the number 
of sites in the fifth dimension. This coupling of 
chiral modes describes a Dirac fermion with mass 
m-ces- In order to have control over the mass of 
the light Dirac fermion, the RH and LH bound- 
aries are explicitly coupled with a weight —am. 
For am ^ am-^es the mixing within the bulk of 
the fifth dimension is negligible. 

The important feature of the domain wall 
(DW) fermion action is that the continuum limit 
a ^ is decoupled from the chiral limit iV^ ^ oo 
at the expense of an extra dimension of size Ns- 
Simulations with DW fermions have shown that 
for lattice spacings roughly 0.1 fm, Ns = 10—20 is 
sufficient for lattice QCD to exhibit continuum- 
like chiral properties, e.g. Ward identities 
and suppression of wrong-chirality operator mix- 
ings @. 

3. RESULTS OF RBC SIMULATIONS 

There are two ingredients which contribute to 
the determination of quark masses from the lat- 
tice: the calculation of hadron masses, decay con- 
stants, and/or matrix elements; and the calcula- 
tion of the quark mass renormalization constant. 

In Refs. 1^,^ the pseudoscalar meson, vector 
meson, and nucleon masses are computed as func- 
tions of the bare quark mass, am. The am corre- 
sponding to the light quark mass is determined by 
extrapolating (aMps/oMy)^ to {M-^jMpf-, and 
then the lattice spacing a can be defined by set- 
ting the p or iV to its physical value. Finally the 
bare am corresponding to the strange quark mass 
is fixed by interpolating the pseudoscalar mass to 
Mk or the vector mass to Mk* or M^. 




am 

Figure 2. Pseudoscalar meson mass squared vs. 
am. 
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Figure 3. Vector meson and nucleon masses vs. 
am. 



Figs. H and H show these hadron masses as func- 
tions of am. Due to a combination of finite spatial 
volume effects, finite effects, and quenching 
errors, the chiral limit is at am = —0.0038(6). 
Interpolating in the hadron masses as described 
above, one finds (with statistical errors only) S: 



bare 
mass 


l/a = 1.91(0.04) GeV 
lattice units MeV 


mi 
ms{<p) 


0.00166(0.00005) 3.17(0.11) 
0.042(0.003) 80(6) 
0.053(0.004) 101(7) 



The scale above has been defined by setting 
the p mass to its experimental value; however 
Mn /Mp = 1.37(5) compared to the experimental 
ratio 1.22, so there is at least a 10% uncertainty 
due to the ambiguity in whether the nucleon or p 
is used to set the scale. This is common among 
quenched lattice simulations with lattice spacings 
C'(O.lfm). 
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In Ref. (2j the renormalization constant is com- 
puted nonperturbatively using the regularization 
independent momentum subtraction (RI/MOM) 
scheme as advocated by the Rome-Southampton 
group (lO|, apphed for the first time to DW 
fermions g. The result at 2 GeV for = 
is 1.63 ± 0.07 (stat.) ±0.09 (sys.) which 
should be compared to 1.32, the one-loop per- 
turbative Z„i | [Tl| , [l^ for these simulation param- 
eters. One inherent assumption of this method is 
that perturbation theory is reliable at the scale of 
the inverse lattice spacing, roughly 2 GeV here. 
As Chetyrkin pointed out |jl^ the convergence 
of the perturbation series for the matching be- 
tween RI/MOM and MS schemes is slow at 2 
GeV. The convergence would be much improved 
for 1/a 2.5 - 3 GeV. 

The result from Ref. |^ at this lattice spacing 
(corresponding to 6/gQ = 6.0) is 



mf ^(2 GeV) = 130 ± 11 ± 18 MeV, 



(1) 



where the first error is the statistical uncertainty 
and the second is the systematic uncertainty. 
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Figure 4. Strange quark mass vs. lattice spacing. 
Except for = point, all iris values computed 
using the p to set the scale, the K to set rUs, and 
the RI/MOM scheme. The — asterisk is an 
average of three recent continuum limit calcula- 
tions. See text for further description. 



4. COMPARISON OF RESULTS 

Fig. ^ shows several quenched lattice calcula- 
tions of TOs plotted vs. lattice spacing squared. 
Except for the point at — 0, the quark masses 
were determined using the procedure described 
above and are plotted with statistical error bars 
only. The square corresponds to the DW fermion 
simulation of Ref. |]. The crosses are from KS 
simulations |]l^, the big (small) diamonds are 
from Wilson (tree-level improved Wilson) simu- 
lations 1 15 1, and the octagon from a nonpertur- 
batively improved Wilson action |l^. Finally, 
the asterisk displays the author's average of three 
recent continuum extrapolations of rUs, as de- 
scribed in the following paragraph. 

Recently there have been three quenched lat- 
tice calculations of nis which use nonpertur- 
bative renormalization and take the continuum 
limit: one with KS fermions ||l^ and two with 
nonperturbatively improved Wilson (NP IW) 
fermions 18|, (see table below). In order to 
precisely compare among different lattice simula- 
tions which contain similar, if not identical, sys- 



tematic effects, it is common practice to quote 
only statistical errors. However, before using the 
results phenomenologically, one must do one's 
best to include these systematic errors. As these 
simulations were performed on large volumes, the 
source of the largest systematic error (within 
the quenched approximation) is the ambiguity of 
which physical quantity is used to set the lattice 
spacing. Ref. Q uses the lattice spacing set by 
the p mass, a(Mp), while Refs. 0,|l^ use a{fK) 
(obtaining identical results with the lattice spac- 
ing set by the hadronic radius parameter rg) and 
report that using a(Mjv) instead would increase 
TTis by 10 MeV. It is commonly seen in quenched 
lattice QCD that 



a(/if) < a{Mp) < a{MN). 



(2) 



Since the latter two references quote a one-sided 
systematic uncertainty, the central value is shifted 
up by half of that uncertainty and the new error 
contains the statistical and systematic uncertain- 
ties added in quadrature. A 5% scale uncertainty 
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is added to the statistical error of Ref. . The 
weighted average of the three results is 106 ± 4 
MeV. 



Ref. 



114] 

1171 
|1S1 



quoted TOs a ^ syst. 



106(7) MeV ±5 MeV 
97(4) MeV +10 MeV 
105(4) MeV +10 MeV 



adjusted 



106(9) MeV 
102(6) MeV 
110(6) MeV 



weighted average 



106(4) MeV 



Of course in nature there are sea quarks, so 
the quenched approximation could prove un- 
reliable. Last year the CP-PACS collabora- 
tion reported deviations from experiment of the 
quenched hadron spectrum in the continuum 
limit at the level of 10% one might expect 
the same level of disagreement for the strange 
quark mass. However, this year the CP-PACS 
collaboration has a preliminary result from two- 
flavor dynamical Wilson fermion simulations of 
mf^(2 GeV) = 84 + 7 MeV pOl. Although 



the mass renormalization has been computed 
only in (improved) perturbation theory, such a 
small strange quark mass is in conflict with lower 
bounds based on the positivity of hadronic spec- 
tral functions [ pT[ |. 

5. CONCLUSIONS 

The RBC collaboration [Q has completed the 
calculation of the strange quark mass at a sin- 
gle gauge couphng G/qq = 6.0 using the do- 
main wall fermion action. The mass renormal- 
ization is determined nonperturbatively through 
the RI/MOM scheme. The result, 130 +11 + 18 
MeV, is consistent others at the same gauge cou- 
pling. Note this is a result at only one lattice 
spacing. Only in the continuum limit can a reli- 
able comparison can be made to the NPSW and 
KS average of 106 + 4 MeV. A final answer for nis 
awaits several unquenched simulations exploring 
all the different systematic uncertainties already 
probed within the quenched approximation. 
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